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We consider estimating the parameter associated with the qubit depolarizing channel when the 
available initial states that might be employed are mixed. We use quantum Fisher information 
as a measure of the accuracy of estimation to compare protocols which use collections of qubits 
in product states to one in which the qubits are in a correlated state. We show that, for certain 
parameter values and initial states, the correlated state protocol can yield a greater accuracy per 
channel invocation than the product state protocols. We show that, for some parameters and initial 
states, using more than two qubits and channel invocations is advantageous. These results stand in 
contrast to the known optimal case that uses pure initial states and a single channel invocation on 
a pair of entangled qubits. 
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I. INTRODUCTION 

Quantum parameter estimation, originally motivated 
by optical metrology [1-3], considers how to estimate pa¬ 
rameters associated with the evolution of quantum sys¬ 
tems by physical means. A theory of quantum estima¬ 
tion, derived from the laws of quantum physics and statis¬ 
tics, has been developed and applied to various evolution 
processes yielding important results for various metrol¬ 
ogy situations [1, 2, 4-11], 

We consider the single qubit depolarizing channel, 
which maps a qubit, initially in the state p , via 

p r & ] p f {\) :=^-^Tr[p\I + \p (1) 

where 0 ^ A ^ 1. The aim of a physical estimation pro¬ 
tocol would be to determine A as accurately as possible 
by subjecting a collection of qubits to this evolution. 

The depolarizing channel is interesting for various rea¬ 
sons. First, depolarization is a standard model of certain 
noise processes and is of general interest for quantum in¬ 
formation processing [10, 12-15]. Specific examples ap¬ 
pear in nuclear magnetic resonance (NMR) [16, 17] and 
optical quantum information processing [18-20]. 

Second, quantum parameter estimation reveals fun¬ 
damental and quantifiable differences with classical ap¬ 
proaches when manifestly quantum resources such as en¬ 
tanglement are used. For example, in the estimation of 
the phase associated with a unitary parameter, classical 
approaches yield an uncertainty that scales as 0(1/^/n) 
where n is the number of times the channel is invoked, 
whereas using entangled states can result in a lower un¬ 
certainty that scales as 0(l/n), thus pointing to a mea¬ 
surable gain provided by entanglement [8, 21, 22], 

Depolarizing channel parameter estimation has been 
investigated previously in the context of finding optimal 
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estimation protocols [23-26] where any possible initial 
state is assumed to be available. It emerges that the op¬ 
timal estimation protocol requires qubit pairs, each ini¬ 
tially in an entangled pure state, and that states display¬ 
ing entanglement or correlations amongst more than two 
qubits will not give any advantages to estimation accu¬ 
racy [25]. 

We consider various estimation protocols for the con¬ 
strained situation where when the available initial state 
of each qubit is mixed. This situation arises in room tem¬ 
perature solution state NMR, where the initial state for 
any nuclear spin is highly mixed [14]. We ask whether, 
for mixed initial state situations, protocols involving cor¬ 
related states can yield gains in estimation accuracy com¬ 
pared to protocols that use product states and whether 
increasing the number of correlated qubits beyond two is 
advantageous. While we will consider these questions 
for situations involving initial states with all possible 
ranges of mixedness, we will focus on the special case 
where the initial state of each qubit is highly mixed. Our 
approach is motivated by similar work for other types 
of channels where gains were found when using corre¬ 
lated states [27, 28] and is also inspired by various other 
approaches to investigating parameter estimation with 
mixed initial states [29-33]. 

This article is organized as follows. Section II offers a 
review of quantum estimation theory. Section III applies 
this to the simplest depolarizing channel parameter chan¬ 
nel estimation protocol, establishes a benchmark against 
which other schemes are to be compared and reviews ex¬ 
isting results for parameter estimation using pure initial 
states. Section IV considers mixed initial state estima¬ 
tion protocols where there are never any correlations be¬ 
tween qubits. Section V, which contains the key general 
results of our article, considers an estimation protocol 
that employs correlated states. Section VI applies these 
to the special case where the available initial states are 
highly mixed. Finally section VII briefly discusses the 
roles of entanglement and quantum discord in the proto¬ 
col that uses correlated states. 
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II. QUANTUM PARAMETER ESTIMATION 

Physical parameter estimation requires that suitably 
prepared quantum systems be subjected to the opera¬ 
tion r(A), which is dependent on a real parameter A that 
is to be estimated. The generic scheme for doing this, 
illustrated in Fig. 1, assumes that m copies of the opera¬ 
tion, or channel, are available. These may be invoked on 
one or more of n systems with invocations surrounded 
by a series of parameter independent unitary transfor¬ 
mations. Each quantum system is prepared in a known 



Pf(A) 


FIG. 1. Quantum estimation procedure using n quantum sys¬ 
tems, each represented by a horizontal line. This procedure 
uses m invocations of the operation r(A), and parameter- 
independent unitary transformations Ui,.Um+i- The ini¬ 
tial state of each quantum system is po and the pre¬ 
measurement state of the collection of systems is pf(A). Mea¬ 
surements on individual systems follow U m +i- 

initial state po- After the final unitary, a measurement is 
chosen and performed on each system and the measure¬ 
ment outcomes are recorded. The final pre-measurement 
state pi (A) may depend on the parameter and thus so will 
the probabilities with which the measurement outcomes 
occur. Thus information about A might be inferred from 
the measurement outcomes by classical data processing, 
thus yielding an estimate for the parameter. 

The fact that measurement outcomes in quantum 
physics are usually statistically distributed implies that 
different runs of the same estimation procedure will yield 
estimates that fluctuate about a mean value. The aim of 
quantum parameter estimation is to tailor the process 
by choosing the initial state po , the parameter indepen¬ 
dent unitaries Ui,..., U m + 1 and the final measurements 
so that the estimates will display minimal fluctuations 
about the actual parameter value. 

The accuracy of the estimation is quantified via these 
statistical fluctuations as follows. Given a set of mea¬ 
surement outcomes, one on each qubit and denoted 
{xi,... ,x n }, an estimate for the parameter A is attained 
via A = A es t(^i, • • •, x n ) where A es t is an estimator func¬ 
tion. Fluctuations in these measurement outcomes from 
one run of the estimation procedure to another gener¬ 
ate fluctuations in A. The extent to which this occurs 
depends on the construction of the estimator function. 
Key requirements for an estimation protocol will be that 
the estimator is unbiased, i.e. (A) = A, and that it mini¬ 
mizes the variance, 


var (A) := ^A — 


( 2 ) 


where the angle brackets indicate the mean over all possi¬ 
ble measurement outcomes. The variance will depend on 
the choice of estimator. However, a key result [9, 34, 35] 
is that the variance on any estimator is bounded from 
below via the Cramer-Rao bound 

™ {X)> ¥(X) (3) 


where 



<91np(:ri,... x„|A)] 2 
dX 


dxi... dx n 


( 4 ) 


is the (classical) Fisher information associated with the 
probability distribution for the measurement outcomes 
p(xi,... x„|A). This is independent of the choice of esti¬ 
mator. Quite generally, there is always an estimator [34], 
which asymptotically saturates the bound of Eq. (3). 
Thus in classical parameter estimation the Fisher infor¬ 
mation quantifies the accuracy of any possible estimate; 
the larger F( A), the better the estimate. 

In classical estimation, it is easily shown [35] that if the 
same estimation trial is repeated m times independently, 
the the Fisher information becomes mF( A), where F( A) 
is the Fisher information for a single trial. Thus the 
Cramer-Rao bound implies a reduction in the variance 
by a factor of m, yielding the typical classical reduction 
in uncertainty of the estimate by a factor of \Jm. 

In quantum evolution processes the probability distri¬ 
bution for the measurement outcomes depends on the 
type of measurement. Thus, given many possible choices 
of measurement, there will be many possible values for 
the classical Fisher information for the same evolution 
process. Quantum estimation aims for a measurement 
which results in the largest classical Fisher informa¬ 
tion. This is constrained and simplified by the quantum 
Cramer-Rao bound, which states that, for any measure¬ 
ment the resulting classical Fisher information satisfies 


F( A) < H{ A) (5) 

where U(A) is the quantum Fisher information and is 
independent of the choice of measurement [5, 9, 35, 36]. 
The quantum Fisher information depends only on pt{\) 
and is calculated via [9] 


H{ A) = Tr /3 f (A)L 2 (A) 


( 6 ) 


where the symmetric logarithmic derivative (SLD) or 
score operator, L( A) is defined implicitly via 


dpf(A) 

d\ 


L(X)pi(X)+pi(X)L(X) 


( 7 ) 


In general the SLD is not immediately apparent from 
Pf(A); sometimes it can be computed using relatively 
simple algebra [28]. However, the SLD and quantum 
Fisher information can always be computed whenever an 
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eigenvalue and eigenstate decomposition for /5f(A) can be 
found. Specifically, if pf(A) = '^jPj\4 > j){ c l ) j\ th en [9] 




j,k 


Pj dX 

(Pj -Pk) 2 

Pj + Pk 


(<Aj 


d\fa) 


dX 


( 8 ) 


Quite generally, there exists a measurement procedure 
which asymptotically attains the quantum Cramer-Rao 
bound [37]. Thus the quantum Fisher information quan¬ 
tifies the accuracy of the quantum estimation process. 
The main task in quantum parameter estimation is to 
provide a protocol giving a final pre-measurement state 
pt(X) with the maximum quantum Fisher information; 
this approach has been widely adopted elsewhere [5, 7 
9, 11, 22, 23, 25-28, 32, 35, 38-42], 

Before continuing, we remark that a useful computa¬ 
tional alternative to Eq. (6) is [28] 


H( A) = Tr 


dpf(X) 


dX 


L(X) 


(9) 


The quantum Fisher information will depend on the 
resources, such as the numbers of qubits, channel invoca¬ 
tions, single qubit unitaries or entangling unitaries that 
are used and various competing estimation protocols can 
only be compared in terms of their resource use. In this 
article we take a common approach, which regards the 
number of channel invocations m in the circuit of Fig. 1 
as the only important resource cost [9, 11, 21, 22, 43]. In 
this context, the goal of any quantum estimation protocol 
is to produce a final state with the maximum quantum 
Fisher information per channel invocation. 

The initial state for the entire system can be assumed 
to be a product of states for the individual systems. If 
necessary, any of Z7i,..., U m +i can be used to generate 
entangled or otherwise correlated states as the procedure 
unfolds. We can then consider two broad classes of esti¬ 
mation schemes. In the first the unitaries Z7i,..., 
never produce any correlation between the individual sys¬ 
tems and the entire system evolves through a series of 
product states. Various classical approaches to estima¬ 
tion, including repetition and averaging, fall into this 
realm. In the second, some of the unitaries Ui ,..., U m +i 
produce states where there is some form of quantum cor¬ 
relation, such as entanglement, between the individual 
systems. Frequently quantum estimation aims to show 
that using such correlated states results in a parameter 
estimate with a greater quantum Fisher information per 
channel invocation [8, 23, 27, 28, 39] than is possible with 
non-correlated states; indirectly this indicates the utility 
of quantum resources for information processing. 


A. Initial state considerations 

Within the framework of Fig. 1, the typical task is to 
find the initial state po and the unitaries Ui,..., U m +i 
that produce a final state /5f(A) that yields the maxi¬ 
mum quantum Fisher information per channel invoca¬ 
tion. The convexity of the quantum Fisher information 
guarantees that this will be attained by using a pure ini¬ 
tial state [23, 40], and most investigations of depolar¬ 
izing channel parameter estimation have considered this 
case [23-26]. However, it may be that case that the avail¬ 
able initial states are not pure. Our aim is to consider 
various estimation protocols when the available initial 
states are limited to those that are mixed. 

The general initial state of a single qubit is 

Po = \ (i + r-«r) (10) 

where r • <r = r x a x + r y a y + r z a z and r x , r y ,r z are the 
components of a three dimensional real vector, r. The 

magnitude of this vector r := \Jr x + r y + r\ is called the 

polarization or purity of the state and it satisfies 0 ^ r ^ 
1. When r = 1 the state is pure and as r, which we use to 
quantify the mixedness or purity of the state, decreases 
the state becomes more mixed. 

Thus the principal issue of such mixed state quantum 
estimation is: given that the available qubits have initial 
states each with polarization r, how do the accuracies, 
quantified by the quantum Fisher information, of various 
parameter estimation protocols compare? Note that, in 
this approach, the orientation of the polarization vector 
r is irrelevant, as this can be modified at no cost with a 
single qubit unitary without changing the polarization r. 


III. DEPOLARIZING CHANNEL PARAMETER 
ESTIMATION: BACKGROUND AND 
BASELINES 


The simplest parameter estimation scheme, which also 
forms the baseline against which all other schemes may 
be compared, is that where the channel is invoked once 
on a single qubit and is described in Fig. 2. This will be 
termed the single qubit, single channel (SQSC) protocol. 



f(A) 


FIG. 2. SQSC protocol. The qubit is initially in initial state 
po. A parameter independent unitary transformation Ui pre¬ 
cedes the single channel invocation, producing a channel input 
state, pi. A measurement follows the channel invocation. 


The unitary Ui of Fig. 2 only changes the orientation 
of r, leaving the polarization unaltered and thus can be 
omitted. Then Eq. (1) yields that the state after one 
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channel invocation is pf = yl + Ar • <rj /2. The quan¬ 
tum Fisher information can be computed by first de¬ 
termining the SLD using Prop. 1 of [28] and then sub¬ 
stituting into Eq. (9). This proposition (appropriately 
corrected) states that if a := Tr(p 2 ) — (Trpt) 2 = 0 


then L = 
then L = 


p dpt _ 3 In |Trp f | -■ 

" a\ a \ pi 


' d\ 


d\ 


p d£f. — 91n H n - 
^ d\ d\ ef 


/ TV pf while if a ^ 0 

/ Tr Pf + ^ln(|a/Trp f |)J. 
Straightforward calculation gives that Tr pf = 1 and 
a = (A 2 r 2 — l)/2. The only circumstances under which 
a = 0 are when r = 1 and A = 1. We henceforth only 
consider A < 1. It follows that, for a single channel invo¬ 
cation on the state of Eq. (10) the SLD is 


L 


A r 2 ~ 1 

1 - A 2 r 2 + 1-A 2 r 2 r ' 


( 11 ) 


Substituting into Eq. (9) then gives that the quantum 
Fisher information for the SQSC protocol using an initial 
state with polarization r is 

= rr^3- (12) 

This is the baseline quantum Fisher information against 
which all other estimation schemes are to be compared. 


A. Pure Initial State Estimation Schemes 


We briefly digress to discuss the previously investi¬ 
gated [23-26] extreme case, where the initial states are 
pure. If it were possible to modify the initial state po¬ 
larization, then Eq. (12), shows that, in accordance with 
general results [23, 40], the largest quantum Fisher in¬ 
formation would be attained for a pure input state, i.e. 
when r = 1, giving 


H, 


s pure 


(A) 


1 

I^A 2 ' 


(13) 


Usually the key issue of interest in quantum parameter 
estimation is whether a greater quantum Fisher informa¬ 
tion can be attained when entangled multiple qubits are 
available as input for the channel. For the depolarizing 
channel it has been shown that entanglement can yield 
a greater quantum Fisher information. Specifically if the 
channel is applied to one of two qubits then the optimal 
quantum Fisher information is attained [23, 26] when the 
two qubits are prepared in a maximally entangled pure 
state, such as |\E'j) = (|00) + 111)) /y/2. The resulting 
quantum Fisher information is 


H, 


pure opt 


(A) = 


(3 + A)(l - A) 


(14) 


The remaining question is whether any additional 
qubits and channel invocations can yield a greater quan¬ 
tum Fisher information per channel invocation. Remark¬ 
ably the result is [25] that the optimal depolarizing chan¬ 
nel parameter estimation scheme uses pairs of qubits, 
each in a maximally entangled pure input state. The 
channel is applied once to only one of each entangled 
pair and the resulting quantum Fisher information per 
channel use is that of Eq. (14). Thus when qubits are 
available in pure initial states then entangling more than 
two qubits cannot yield any advantages for depolarizing 
channel parameter estimation. However, we stress that 
this does not necessarily pertain to the situation that we 
address, in which the available input states are not pure. 


IV. MULTIPLE CHANNEL USES: 

UNCORRELATED STATE PROTOCOLS 

The previous discussion motivates the central ques¬ 
tions of our work: given a collection of qubits, each ini¬ 
tially in a state with the same polarization r, is it possible 
to enhance the parameter estimation accuracy by using 
multiple channel invocations or entanglement or correla¬ 
tions amongst more than two qubits? If so to what extent 
can the accuracy, as measured by the quantum Fisher in¬ 
formation per channel use, be enhanced? 

A first strategy for enhancing estimation accuracy is 
one where the same process is repeated independently, as 
is done for classical repeat and average schemes. Specif¬ 
ically we consider the protocol of Fig. 3, called the inde¬ 
pendent channel use protocol. Here there are m qubits, 
each in the initial state po and subject to a parameter in¬ 
dependent single qubit preparatory unitary. The channel 
is subsequently invoked exactly once on each qubit. In 
this protocol there are never any correlations between the 
qubits during the entire process; in this sense it resembles 
a classical approach. 



A (A) 


FIG. 3. Independent channel use protocol. There are m single 
qubits, each initially in the state po. Qubit j is subjected 
to a single qubit unitary, Uj, which is independent of the 
channel parameter, producing the input state pi. Thereafter 
a depolarizing channel is applied once to each qubit. 

The quantum Fisher information is additive in such 
cases [28, 35] and thus the quantum Fisher information 
for the independent channel use protocol with m channel 
invocations is 


mr 2 


It is straightforward to show that H pure opt (A) A 
H s pure (A) whenever A > 0. 


H( A) = mH s { A) 


1 — A 2 r 2 


(15) 


















5 


Although this is an improvement over the SQSC proto¬ 
col, clearly the quantum Fisher information per channel 
use in the independent channel use protocol is identi¬ 
cal to that of the SQSC protocol. Thus the independent 
channel use protocol attains the same estimation accu¬ 
racy, measured in terms of quantum Fisher information 
per channel invocation, as the SQSC protocol. 

A second strategy, illustrated in Fig. 4, is that where 
only one qubit is available and the channel is invoked 
sequentially a total of m times. Successive channel in¬ 
vocations are interspersed with parameter-independent 
unitaries. Assuming the standard initial state, po, the 


4 

3 

Gseq(A) 2 
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FIG. 5. Sequential channel use gain for m = 3. 
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FIG. 4. Sequential channel use protocol. A single qubit, ini¬ 
tially in the state po is subjected to m channel invocations. 
Each channel invocation is preceded a single qubit unitary, 
U\ ... Um, all of which are independent of the channel param¬ 
eter. 


effect of U\ is to rotate the polarization vector, r —!> iq 
without altering its magnitude. Then Eq. (1) implies 
that the state after the first channel invocation is pi = 
(/ + Ari • <r)/2. Using the same argument, the state after 
the second channel invocation is p 2 = (I + A 2 r 2 • <r)/2 
where U -2 rotates iq —> r 2 - Continuing, the state prior 
to measurement is pf = (I + A m rf ■ &)/2 where ry = r. 
A calculation similar to that giving Eq. (12) yields that 
the quantum Fisher information for m sequential channel 
invocations on a single qubit is 

^e q (A)=m 2 ^^. (16) 

Thus the quantum Fisher information per channel invo¬ 
cation is 

y2m—‘2 r 2 

-^seq/channel(A) = m ^ _ ^2m f 2 ' 


To compare this to the SQSC protocol we define the se¬ 
quential channel use gain as 


Gse q (A) := 


-^seq/channel ( A) 


^ _ ^2rrij.2 


(18) 


Figure 5 illustrates the typical gain, plotted here for the 
m = 3 case. This shows that the sequential channel use 
protocol is advantageous for certain ranges of parameter 
values and polarizations but disadvantageous for others. 


The sequential channel use gain decreases monotoni- 
cally with r and increases monotonically with A (these 
are proved in Appendix A). A more detailed analysis of 
this separates the pure initial state case, i.e. r = 1, from 
mixed initial state cases, i.e. r < 1. For the pure initial 
state case, it is shown in Appendix A that 

0 < G S eq(A) < 1. (19) 


For the mixed initial state case, it is also shown Ap¬ 
pendix A that 

0 sj G S eq(A) ^ m. (20) 

In both cases the upper bound is attained when A = 1. 

It follows that if the available initial state is pure there 
is no advantage to the sequential channel use protocol. If 
the available initial state is mixed then, for certain initial 
polarizations and the parameter values, there can be an 
advantage to the sequential channel use protocol. 

It remains to assess whether for mixed initial states, 
there is any advantage to increasing the number of chan¬ 
nel invocations. To do so, we compare the quantum 
Fisher information per channel use of Eq. (17) for m + 1 
channel invocations versus that for m channel invoca¬ 
tions. An advantage is attained when the ratio of that 
for m + 1 to that for m exceeds 1, and straightforward 
algebra shows that this is equivalent to r 2 ^ [A 2 (m + 1) — 
to]/ X 2m+2 . The right hand side of this increases mono¬ 
tonically with A and attains a maximum of 1 when A = 1. 
Also it is negative when A < y/m/(m + 1), implying that 
when A < s/m/(m + 1) there is no polarization for which 
an additional channel use will give an advantage. Thus 
for any polarization there will be some parameter values 
for which additional invocations in the sequential channel 
use protocol will be advantageous. However, the range of 
possible parameter values f or which thi s occurs, bounded 
from below by at least s/m/ (to + 1), shrinks with in¬ 
creasing channel invocations. 


V. MULTIPLE QUBIT OR CHANNEL USES: 

CORRELATED STATE PROTOCOL 

An additional resource available to multiple qubit sys¬ 
tems, is the existence of states in which there are correla¬ 
tions between the individual qubits and we now consider 
whether these can enhance the accuracy of parameter 
estimation for the depolarizing channel. 

We consider the particular correlated state protocol 
procedure of Fig. 6; this was previously investigated for 
other channels [27, 28]. We assume that each qubit is 
initially in the particular state po = (/ + ra y )/2 with 
polarization r and that a preparatory unitary U prep , il¬ 
lustrated in Fig. 7, is subsequently applied jointly to the 
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qubits. Denote the density operator for the system af- 





> PfW 
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FIG. 6. Correlated state protocol considered in this article. 
There are n single qubits, each initially in the state po. These 
are subjected to a joint preparatory unitary, U prep , producing 
the channel input state pi. Thereafter a depolarizing chan¬ 
nel is applied once to each of the uppermost m qubits. The 
preparatory unitary is illustrated in Fig. 7. 


ter the preparatory gate by p\ := U piep piU2 ep . It can be 


The diagonal terms are 

dj = 2^+1 [(1 + rf(l-r)^ 

+ (1 + r) n ~ j {1 - r) j 


(24) 


As described previously [28], this scheme results in a den¬ 
sity operator which has supports consisting of mutually 
orthogonal two dimensional subspaces; in terms of ma¬ 
trix representations it can be expressed as a sum of 2 x 2 
matrices acting on distinct vector spaces. 

When the initial states are pure, i.e. r = 1, the di¬ 
agonal and counter-diagonal terms are all zero except 
Co = — c n = —1/2 and do = d n = 1/2 and pi is the 
density operator corresponding to the GHZ-type state 
(|0...0) + |l...l))/v / 2. 

Suppose that the channel acts once on each of m qubits 
corresponding to the least significant (or rightmost) dig¬ 
its in the computational basis representation. The action 
of the channel on multiple qubits is readily described by 
noting that Eq. (1) implies that, for a single qubit, 



FIG. 7. The preparatory gate, U plep , consists of a collection of 
controlled-!? operations applied to each distinct pair of qubits 
followed by a Hadamard transformation applied to each qubit. 

shown [28] that 

(JV—1)/2 

Pi = y ' Px ( 21 ) 

x=0 


\ x )(v\ ^ A|a;)(y| and (25a) 

|aj)(a;| p\x){x\+q\x)(x\ (25b) 

where |x) and | y) are computational basis states such 
that x y and x := 1 — x. Here p := (1 + A)/2 and 
q := (1 — A)/2. A general state of multiple qubits can 
be expressed in terms of linear combinations and ten¬ 
sor products of basic operators such as \x){y\ and thus 
Eqs. (25) can be used to describe the action of any num¬ 
ber of channels acting on multiple qubits. 

Eq (25a) implies the counter-diagonal terms of the ma¬ 
trix representing pi (A) are attained from those of pi by 
multiplication by A m . Eq. (25b) implies that the diagonal 
terms of the matrix representing pt(X) result from rear¬ 
rangements of those of pi multiplied by various factors 
of p and q. The remaining terms are again all 0. An 
additional simplification, demonstrated in appendix B is 
that (N — x\pi(X)\N — x) = {x\pi (A)|a;). Thus the final 
density operator has a form similar to that for pi, i.e. 


where N = 2 n — 1 and 

Px =dj(x) l a; )( ;l: |+|A’ — x)(N — x\ 

+ icj( x ) \x)(N - x\-\N - x)(x\ 


( 22 ) 


The counter-diagonal terms of the resulting matrix of 
Eq. (22) are constructed by first computing the number 
of zeroes j(x) in the binary representation x n ... X\ for x 
and then using 


1 


3 2n+l 


d = —- [(1 + r) j (1 - r) n ~ j 

- (1 + r) n ~ j (l - r) j 


(23) 


(JV-l)/2 

Pf(A) = ^2 P {x 

x=0 


(26) 


where 


Pf x —df 


|cc)(ar|-F|A r — x)(N — x\ 


■ ict x |a:)(./V — x\—\N — x){x\ 


(27) 


Eq. (25a) yields 


Cfx = A m Cj( x ) ( 28 ) 

for the counter-diagonal terms where j (x) is the number 
of zeroes in the bit string for x. 
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Determining the diagonal terms is further simplified 
by noting that the preparatory unitary is invariant un¬ 
der interchange of any two qubits and that this is also 
true of the channel operation provided that the channel 
either acts on each qubit that is interchanged or that it 
acts on neither of the qubits that are interchanged. The 
invariance of the preparatory unitary is readily apparent 
from its structure, as shown in Fig. 7. More precisely, 
suppose that Syj is the unitary operation which swaps 
the qubits labeled j and j'. Then if 1 ^ j < j' ^ to 
( recall that the channel acts on the rightmost m bits), 

(y\S],jPiS rj \x) = (y\f>f\x) (29) 

with an identical result when to + 1 Sj j < j' ^ n. 
Thus if the bit string for x is x n ... Xm+iXm ... xi, then 
d{ x = {x\pi\x) is invariant under the interchange of 
any two of {x n ,..., x m+ i} and is also invariant under 
the interchange any two of {x m ,...,xi\ (it is not in 
general invariant under interchanges between one cle¬ 
ment from one of these sets and another from the other 
set). For example, if the channel is applied the two 
least significant qubits out of a total of five qubits then 
df 6 = df 5 (or, using a binary representation for the in¬ 
dex subscript, dfoono = dfooioi) and d{ 2 6 = df 2 i (or 
df noio = dfioioi)- A consequence of this is that df x 
only depends on the number of zeroes in the leftmost 
n — m bits of the binary representation of x and also the 
number of zeroes in the rightmost m bits of the binary 
representation of x. This consideration and a derivation, 
provided in appendix C, yields 

m ^max / \ / \ 

dfx=E^ m_fc E (i)uin ( 3 °) 

k =0 ' / V / 


where u = u{x) is the number of zeroes in the leftmost 
n — to bits (i.e. x n ... x m +i ) of the binary representation 
of x, v = v(x) is the number of zeroes in the rightmost 
to bits (i.e. x m .. .Xf) of the binary representation of £, 
/ m in = max (k + v — to, 0), and Z max = min (k, v). Clearly 
df x depends on u, v and m and A, via p = (1 + A)/2 and 
q = (1 — A)/2. Since A is implicit throughout this work 
and to is fixed by the number of channel invocations we 
will find it convenient to express df x as df x = df (u,v). 

According to Eqs. (26) and (27), pf is a sum of op¬ 
erators, each acting on a two dimensional space, which 
have mutually orthogonal supports. Thus Prop. 2 of [28] 
implies that the quantum Fisher information is 


this and Eq. ( 8 ) yields 


77 C ott (A) 


(AT—1)/2 

= E 


x =0 


dj 


2 


A 2m c 2 f s 

J( x ) 


x [dfx (d f ^ + m 2 A 2m 2 c 2 (a:) ) 
—2rnA 2 m_ 1 d f x£j( x ) 


(32) 


where the dot indicates differentiation with respect to 
A. Although the sum in Eq. (32) is over x, the terms 
indexed as such depend on u and v and it follows that 
the quantum Fisher information can be expressed as a 
sum over these. In appendix D it is shown that if m < n 
then 


m n—m 


idcorr(A) =EE 


v=0 u =1 


n — m — 1 
u — 1 


rf2 _ \2m r 2 

A c u+v 


j ( j2 | *2 \ 2 m—2 2 \ 

+m A c u + v J 


— 2 toA 2 


1 df c: 


u-\-v 


(33) 


where df = df (w, v) is given by Eq. (30). If n = m then 


v—1 ^ 


1 / d 2 — A 2m cl 


[df [d: 


d? + to 2 A 2m 2 c 2 


—2m\ 2m ~ 1 df cl 


(34) 


where df = df (0,w) is again given by Eq. (30). In both 
cases Cj is given by Eq. (23). These represent the sim¬ 
plest expressions that we have been able to obtain for 
the quantum Fisher information for the correlated state 
protocol. 

The quantum Fisher information per channel invoca¬ 
tion is 


gcorr/channel(A) = ^ ^ ■ (35) 

TO 

We aim to compare the performance of the correlated 
state protocol to that of the SQSC protocol in terms 
of the quantum Fisher information per channel use and 
given that the initial qubit polarizations are the same. 
We define the correlated state protocol quantum Fisher 
information gain as 


£? C orr(A) : — 


hi car r / channel ( A ) 

hJ\) 


(36) 


(JV—1)/2 

-f^corr(^) = E H x { A) (31) 

31=0 

where H x ( A) is the quantum Fisher information associ¬ 
ated with pf x . A derivation (see appendix D) based on 


and whenever this exceeds 1 , the correlated state protocol 
outperforms the SQSC protocol. 

The expression for the quantum Fisher information is 
a rational polynomial in r and A and even with as few as 
two qubits, this is too complicated to yield insight into 
the estimation accuracy for the correlated state protocol 
for all values of polarization and the parameter. While 











Gain for n = 4, m — 4 


the same applies to the gain, aspects of its behavior can 
be illustrated graphically. 

First, suppose that there is only one channel invoca¬ 
tion (m = 1). Fig. 8 illustrates the gain for the examples 
n = 2 and n = 5. These suggest that when the channel 


Gain for n = 5, m — 1 



Gain for n = 2, m = 1 



FIG. 8. Gain using one channel invocation. Both plots are 
over the range 0.05 ^ A ^ 0.95 and 0 ^ r ^ 1. 


is invoked once only, the correlated state protocol will 
be more accurate than the SQSC protocol, assuming the 
same initial state polarization for both. However, a gen¬ 
eral proof of this is elusive; the methods of [28] used to 
prove a similar fact for the phase-flip channel are not 
immediately applicable here. 

Second, suppose that there are multiple channel invo¬ 
cations (to > 1). Fig. 9 illustrates the correlated state 
gain of Eq. (36) for various numbers of channel invoca¬ 
tions on 4 qubits. These show that when the channel 
is invoked multiple times, there are polarizations and pa¬ 
rameter values for which the the correlated state protocol 
will be more accurate than the SQSC protocol but there 
are also parameter values for which it is inferior. 

Given that the channel is invoked multiple times, the 
correlated state protocol can also be compared to the se¬ 
quential channel protocol that uses the same number of 
channel invocations. We define the correlated vs sequen¬ 
tial channel gain as 


Gcorr/seq(A) ■ 


-ffcorr/channel (A) 
-^seq/channel (A) 


(37) 


where the number of channel invocations is the same for 
the two protocols. Figure 10 suggests that when the 
channel is invoked a fixed number of times, the corre¬ 
lated state protocol is more accurate than the sequential 



Gain for n = 4, m = 3 



Gain for n = 4, m = 2 



FIG. 9. Gain using multiple channel channel invocations on 
4 qubits. All plots are over the range 0 ^ A, r < 0.95. 


protocol, using the same number of channel invocations, 
over all parameter values and for all polarizations. 


VI. WEAK INITIAL POLARIZATION CASES 

In some situations the initial state polarization is very 
small, i.e. r€l. This is true of room-temperature solu¬ 
tion state NMR [14], where typically r « 10 -4 . 

When r<l the SQSC protocol quantum Fisher in¬ 
formation becomes 

Hs(A) « r 2 (38) 

and the sequential channel use protocol quantum Fisher 
information per channel use becomes 

iWhannel(A) ~ mX^^r 2 , (39) 

each to lowest order in r. In appendix E we show that the 
correlated state quantum Fisher information per channel 
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Gain for n = 4, m = 3 Gain for n = 4, m = 3 



FIG. 10. Correlated vs sequential channel gain using multiple 
channel channel invocations on 4 qubits. Plots on the left are 
over the range 0.05 ^ A, r ^ 0.95. Plots on the right provide 
details for A ^ 0.5. 


use becomes 


#corr/channel(A) « mn\ 2m V, (40) 

to lowest order in r. 

Here the sequential channel use gain approximates to 

G seq (A) « m\ 2m ~ 2 . (41) 

Eq. (41) shows that the gain exceeds 1 when the depo¬ 
larizing channel parameter satisfies A ^ A cu toff(m) where 
A cu toff(m) := m 1 /( 2_2rn ) is the characteristic parameter 
cutoff. The characteristic parameter cutoff, plotted in 
Fig. 11, is a monotonically increasing function of m, with 
a minimum value of e -1//2 = 0.61, attained as m —> 1 and 
a maximum value of 1 attained as m —> oo. This indi¬ 
cates that for any given parameter value, there will be a 
maximum number of channel invocations beneath which 
the sequential channel use protocol offers enhanced pa¬ 
rameter estimation accuracy and beyond which it will 
not do so. Figure 11 shows that as A increases, accu¬ 
racy gains will persist with increasing number of channel 
invocations. Note that the larger A is, the weaker the de¬ 
polarization is, and for such weak depolarization cases en¬ 
hanced estimation accuracies can be attained with mul¬ 
tiple sequential channel invocations. Thus, for low polar¬ 
izations, the sequential channel use protocol can yield a 
greater estimation accuracy per channel invocation than 
the SQSC protocol. However, this requires that the pa¬ 
rameter value exceed the characteristic parameter cutoff 
and Fig. 11 indicates that such advantages will be absent 
for low enough parameter values and, even when there is 
an advantage, increasing the number of channel invoca¬ 
tions will eventually eliminate this advantage. 

The correlated state protocol gain (versus the SQSC 
protocol) approximates to 


2m —2 


(42) 



FIG. 11. Parameter cutoff function and its square versus m. 


We consider two cases separately. 

First, with only one channel invocation the correlated 
state protocol offers an n-fold gain in accuracy regardless 
of the parameter value, for low initial polarization states. 
This is strikingly different from the pure state case where 
the optimal accuracy is attained [25] with a single channel 
invocation on one of only two qubits (n = 2 and m = 1); 
by contrast when the initial state polarization is low and 
the channel is only invoked on one qubit, then adding 
more qubits always increases the estimation accuracy. 

Second, with multiple channel invocations, Eq. (42) 
indicates that enhanced estimation accuracy will occur 
if A ^ (TOn) 1 ^ 2-2 " 1 ) = A cu t 0 ff(m) n 1 /( 2-2m ). The right 
hand side forms a new cutoff which is again an increasing 
function of the number of channel invocations but a de¬ 
creasing function of the total number of qubits. Thus in¬ 
creasing the number of spectator qubits while leaving the 
number of channel invocations fixed extends the range of 
parameter values over which the correlated state proto¬ 
col is advantageous. Furthermore Eq. (42) shows that 
adding more spectator qubits while keeping the number 
of channel invocations fixed will increase the quantum 
Fisher information per channel use; the increase can be 
arbitrarily large by adding sufficiently many spectator 
qubits. Unlike the pure initial state case, for low polar¬ 
ization initial states, there will always be an advantage 
to adding additional spectator qubits. 

Separately, when m < n we can assess whether 
one additional channel invocation yields any advantage. 
Eq. (40) implies that starting with m channel invo¬ 
cations and adding one more invocation increases the 
quantum Fisher information if m < A 2 /(l — A 2 ), de¬ 
creases it if m > A 2 /(1 — A 2 ) and makes no differ¬ 
ence if to = A 2 /(1 — A 2 ). Thus the optimal number of 
channel invocations in the correlated state protocol is 
Too p t = LA 2 /(1 — A 2 )J + 1, unless A 2 /(l — A 2 ) is an inte¬ 
ger, in which case, m op t = A 2 /(l — A 2 ). Examples of this 
and the associated gains are provided in Table I. When A 
is close enough to 1, this offers considerable opportunity 
for improved accuracy. Again unlike the pure initial state 
case, for low initial polarization, increasing the number 
of channel invocations in the correlated state protocol can 
result in enhanced estimation accuracy, per channel use. 

In many physical situations every qubit is subjected to 


Gcorr(A) ~ mnX 
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A 

ITLopt 

Optimal gain 

0.700 

l 

l.OOn 

0.800 

2 

1.28n 

0.900 

5 

2.15n 

0.950 

10 

3.97 n 

0.990 

50 

18.67 n 

0.995 

100 

37.07 n 


and thus when r« 1, the correlated state protocol will 
always offer a better estimation accuracy than the se¬ 
quential channel protocol provided that the number of 
channel invocation for each is the same. 


VII. QUANTUM CORRELATIONS IN THE 
CORRELATED STATE PROTOCOL 


TABLE I. Optimal number of correlated state protocol chan¬ 
nel invocations and optimal gain, i.e. Eq. (42) evaluated at 
m 0p t, whenever m < n. 


the depolarizing channel. When the channel is invoked 
on all qubits of the correlated state protocol (m = n), the 
gain becomes G cor r(A) ~ m 2 A 2m_2 and this implies that 
the correlated state protocol provides enhanced accuracy 
whenever A ^ ^cutoff ( m )- This lower bound is plotted 
in Fig. 11. This is again an increasing function of to, 
approaching e -1 = 0.37 as n = m —> 1. The slightly 
lower cutoff implies that there is somewhat greater scope 
for enhanced parameter estimation when all qubits are 
subjected to the channel. However, for a given parame¬ 
ter value, this cutoff establishes an upper bound on the 
number of qubits/channel invocations for which the cor¬ 
related state protocol provides any gain. 

We can again ask when adding one more qubit and 
channel invocation gives any advantage. Equation (42) 
with m = n implies that if there are originally to channel 
invocations then adding one more gives a greater quan¬ 
tum Fisher information per channel use if m ^ m op t 
where TO op t = [A/(1 — A)J + 1 unless A/(l — A) is an 
integer, in which case, m op t = A/(l — A). Examples of 
this and the associated optimal gains are provided in Ta¬ 
ble II. Clearly, for relatively weak depolarizing channels, 


A 

ITLopt 

Optimal gain 

0.500 

l 

1.00 

0.700 

3 

2.16 

0.900 

9 

15.01 

0.950 

19 

56.96 

0.970 

33 

155.03 

0.990 

100 

1367.00 


TABLE II. Optimal number of correlated state protocol chan¬ 
nel invocations and optimal gain, i.e. Eq. (42) evaluated at 
TOopt, whenever the channel is invoked once on each qubit 
(to = n). 

i.e. A « 1, substantial enhancements of estimation accu¬ 
racy are possible with the correlated state protocol. 

Finally, the correlated state protocol and sequential 
channel protocol can be compared in the low initial po¬ 
larization limit, given the same number of channel invo¬ 
cations. Here the correlated vs sequential channel gain 
as 


Gcorr/seq(A) ~ Tl 


(43) 


The presence of quantum correlations such as entan¬ 
glement or discord in the correlated state protocol may 
be assessed for the n = 2 along the same lines as was 
done for the Pauli channel parameter estimation of [28]. 
Then Eqs. (27), (28) and (30) yield that with to channel 
invocations, 


Pf= 4 


fl + A m r 2 0 0 

0 1 - A m r 2 0 

0 0 1 - A" 

\ —2ir\ m 0 0 


2ir\ r 

0 

0 

1 + A m 


(44) 

in the computational basis. The state of the system prior 
to channel invocation can be obtained from this by set¬ 
ting A = 1. 

The presence of entanglement can be assessed by con¬ 
sidering partial transpose of pi and this is 


d PT - - 
Pl ~4 


/1 + A" 
0 
0 

V 0 


0 

1 - A m r 2 
-2irX m 
0 


0 

2ir\ n 
1 - A m 
0 


0 

0 

0 

1 + A” 


(45) 

The state is separable if /5 PT is positive [44, 45]. The 
eigenvalues of this will only be positive if and only if 
0 ^ r ^ \J\ + 1/A m — 1. Thus the state prior to channel 
invocation will be separable if and only if 0 < r < y/2 — 1. 
The correlated state protocol can give estimation gains 
as r —> 0, and this clearly corresponds to a state which is 
separable throughout the protocol. Thus the correlated 
state protocol gains cannot always be ascribed to the 
presence of entanglement. 

Quantum discord describes correlations between quan¬ 
tum systems in terms of mutual information and mea¬ 
sures the extent to which these exceed any possible clas¬ 
sical correlations between the systems [27, 46]. It is pos¬ 
sible to compute the quantum discord for states with the 
form [47, 48] 


P 


1 

4 


3 =1 


(46) 


where Cj are all real and a\ = a x , ■ • ■ ■ Defining 


Ao 

:= 1 - Ci - c 2 - c 3 

(47a) 

Ai 

:= 1 — Ci + c 2 + c 3 

(47b) 

A 2 

:= 1 + ci - c 2 + c 3 

(47c) 

A 3 

:= 1 + ci + c 2 - c 3 

(47d) 
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and 


c := max{|ci|, |c 2 |, |c 3 |} (48) 


yields that the quantum discord is [47] 


1 ^ 1 c 

Q(p ) = 4 5Z l0 §2 - lo §2 (1 - c) 

1=0 

log 2 (l + c). (49) 

Following the scheme of [28], the final density operator 
can be brought into the form of Eq. (46) by applying the 
single qubit unitary 


U = 


( 0 e i7r / s \ 

^g-ivr/8 0 J 


(50) 


to each qubit and this will not alter the discord. Then 


A 0 = 1 — A m r 2 (51) 

Ai = 1 + 2r\ m + A”Y 2 (52) 

A 2 = 1 - 2rX m + A m r 2 (53) 

A 3 = 1 - A m r 2 (54) 


c = max {A”V 2 , A m r} = A m r (55) 


since 0 < r < 1. Again the discord for the correlated state 
protocol prior to channel invocation can be determined 
by setting A = 1 and this gives 

Q(p) = i°g 2 (1 + 0 + ~2~ lo g 2 (1 — r )- ( 56 ) 


Then ^ > 0 whenever r > 0 and Q = 0 when r = 0 
implies that the state prior to channel invocation has 
non-zero discord whenever the polarization of each qubit 
is non-zero. In this sense, we can regard the correlated 
state protocol as using states that have quantum corre¬ 
lations. 


VIII. DISCUSSION 

We have considered estimation of the depolarizing 
channel parameter for various situations where the avail¬ 
able initial state of each qubit is mixed and we quantified 
the accuracy of each in terms of the quantum Fisher in¬ 
formation after channel invocation. We showed that for 
certain ranges of initial state polarization and channel 
parameter, a protocol that uses correlated states prior to 
channel invocation attains a larger quantum Fisher in¬ 
formation per channel invocation than either a protocol 
which uses product states throughout or another proto¬ 
col which uses a sequence of channel invocations on a 
single qubit. For qubits with small initial polarizations 
we showed that using more than two qubits and more 


than two channel invocations can be advantageous; this 
is in contrast to the situation where qubits in pure initial 
states are available. 

Our results suggest that, as with Pauli channel pa¬ 
rameter estimation, a single invocation of the channel on 
many correlated qubits gives gains for all parameter val¬ 
ues and initial polarizations. We were only able to prove 
this for the depolarizing channel in the limit as the initial 
polarization approaches zero, but we were also unable to 
generate a counterexample. This was not central to our 
results but it does merit further investigation. 

The source of the gain cannot always be attributed to 
the presence of entangled or non-separable states prior 
to channel invocation; we demonstrated this for the two 
qubit example. We also showed that there is non-zero 
discord present prior to channel invocation in this exam¬ 
ple. A full analysis of the evolution of the discord and its 
possible relationship to gains in the quantum Fisher in¬ 
formation might provide insight into the source of these 
accuracy gains; this is beyond the scope of our article 
and remains to be investigated. 

Finally, we chose one particular preparatory unitary 
transformation for the analysis of the correlated state 
protocol; that this yields gains in some cases is interesting 
enough. However it is unclear whether this is optimal 
or whether there exists another parameter independent 
unitary that offers greater gains. 


Appendix A: Sequential channel use gain 

The behavior of G seq (A) with respect to r can be as¬ 


sessed by first setting x := r 2 . Then 


G seq (A) — m 


1 — X 2m x 


gives that 


dGseq(A) 

dx 


= —mA 


2m 


1 + A 


2m—2 


(1 - X 2m x) 


<0. 


(Al) 


(A2) 


Since r increases as x increases this proves the fact that 
Gseq(A) is a monotonically decreasing function of r. 

To assess the behavior of G seq (A) with respect to A, let 
y := 1/A 2 and note that 1 ^ y. Then 


G S eq(A) — TTl 


y-r 

ym _ j,2 


(A3) 


and 


m —1 

„ 2 / 2 \ \ ^ „ k 2m— 2k—2 

y -r ={y-r)) j yr 
fc=o 

+r 2m -r 2 . 


(A4) 


Thus 


G S eq(A) = -,-. 

Y^Qy k r 2m ~ 2k ~ 2 + {r 2m — r 2 )/{y — r 2 ) 


(A5) 
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The term ( r 2m — r 2 )/(y — r 2 ) is negative since 0 ^ r ^ 1 
and 1 ^ y and thus decreases as y increases. The re- 
maining term 2^fc=o y r increases as y increases. 

Thus G seq (A) decreases as y increases. However, y de¬ 
creases as A increases and this proves the fact that 
G S eq(A) is a monotonically increasing function of A. 

For the pure initial state case, r = 1 and Eq. (18) gives 

G se q(A) = m y m _ 1 (AS) 

where y = 1/A 2 . Then 


acts on qubit 1 and then on qubit 2. This maps the 
diagonal component of the prepared density operator as 

f(i) ^ ^ a 

Pprep diag 1 ^ PPprep diag H"~ (?-^lPprep diag-^1 

f (2) 2 - 

1 ^ P Pprep diag 

~\~PQ [^lPprep diag-^1 H - -^2Pprep diag-^ 2 ] 

+q 2 X2Xip plep diag-A' 2^1 (B4) 

where X 2 '■= This extends to action of 

the channel once on each of qubits 1,2,3,.. .to and gives 


gives 


m —1 

y m -l = {y-l)Y j y k 


k —0 


G seq ( A) = 


r™" 1 yk ' 


(A7) 


(A8) 


This is a monotonically decreasing function of y ^ 1 
and thus attains a maximum of 1 at y = 1. Thus 
0 ^ G seq (A) ^ 1 for a pure initial state, where the lower 
bound stems from the fact that the gain is a ration of 
positive quantities. 

For the mixed initial state case, G seq (A) is a monoton¬ 
ically increasing function of A and attains a maximum of 
to when A = 1. This shows that 0 ^ G seq (A) ^ to when 
r < 1. 


Appendix B: Symmetries in the diagonal terms of pf 


Eqs. (25) give that the diagonal terms of pf only arise 
from the action of the channel on 


Pprep diag 


(JV—1)/2 

x—0 


| x) (x\ + 1 N — x) (N — x | 


(Bl) 


and that this will be mapped to the diagonal component 
of the final density operator, denoted pt diag and defined 

via /5 prep diag r0 hA sr p { diag . A useful expression for p t d i ag 
involves the following representation of the channel on 
any diagonal operator. Eq. (25b) gives 


p 

\x)(x\>-^ p\x)(x\+qa x \x)(x\a x . (B2) 

Thus if the channel acts on the rightmost qubit, labeled 1, 
then its action on the diagonal component of the density 
operator is 


Pprep diag 1 ^ Pi diag 


= E«V”“' X 


k —0 


'y ' Xlk ■ ■ ' X-hPpiep diag Xl k ' ' ' Xh (B5) 


where the inner sum is over all (Zfe,..., Zi) that satisfy 
Ik > Ik- 1 > • • • > h and range over the entire possi¬ 
bilities starting with l\ = 1 to Ik = m. If k = 0 the 
inner sum is taken to mean p prep diag- For example, if 
k = 3 and to = 4 there are three indices I 3 > I 2 > h 
and the summation range yields the following possibili¬ 
ties for (Z 3 , Z 2 , Zi): (3, 2,1), (4,2,1), (4, 3,1), (4,3, 2). This 
consists of the set of all flips on three out of the four 
qubits with least significant label. A more symmetrical 
alternative is 


m 1 

k —0 

y ' Xlk ' ' ' X-hPprep diag Xl k ' ' ' Xh (B6) 

where the sum is over all to ^ Ik, Ik- 1 , - ■ •, Zi ^ 1 where 
no two indices are equal. We introduce the following no¬ 
tation to render this more compact. Let l := (Ik, ■ ■ ■ ,li) 
where, to ^ Ik, Ik- 1 , ■ • ■ ,h -S 1 and no two elements in 
this list are equal. Then define Xp := Xi k ■ ■ ■ Xh . It fol¬ 
lows that 


qkpm-k 

Pi diag = Af Pprep diagAf ■ (B7) 

/c=0 ' f 

We use this to prove that 

(N - x\pt(X)\N - x) = (z|/5f(A)|x). (B8) 

To do so, note that | N—x) = X n ■ ■ ■ X\\x) and let A to t := 
X-n ■ ■ ■ Xi- Thus, using X% ot = X tot , 


Ppiep diag ^ PPprep diag + qXlpprep diagAf 1 (B3) ( N ai|Pf(A)|iV x) - (N ^|diagl^V X) 

A A A — (ic| Xf 0 fPf diagAtot |-^) ■ (B9) 

where X\ := I ® • • • 0 I ® a x and the superscript on T 
indicates action on qubit 1. Suppose that the channel Eq. (B7) gives 







13 


TYl Ytl/ _/c 

(N - x\p { (X)\N -x) = <*|^ X tot X r p pTep dias X r X tot \x) 

fc=0 ' i 

m m —k 

= 55 — 55^|Xrl to tPp r ep diag*tot*f |*> (BIO) 

fc=0 ' f 


since X to t commutes with Xj- . Then Eq. (Bl) gives 
AtotPprep diagAtot — Pprep diag- This and Eq. (BIO) give 

(N -x\p { (X)\N-x) = (x\± ^ x 

k =0 

5 ' Xf Pprep diag-Xf I*) 
f 

= <ar|/3 f (A)|a:). (Bll) 

I 


d {x = '£q k p rn ~ k 55 (x n ...x lk ...x h 

k =0 


with the same constraints for the inner sum as for 
Eq. (B5) and where Xi = 1 — xi. We shall evaluate the 
inner sum for any value of k, which represents the num¬ 
ber of flips on the bits within the binary representation 
of x. 

The fact that the diagonal terms of pf only depend on 
the number of (and not their order) zeroes, u, in the n—m 
most significant bits and also the number of zeroes, v, in 
the m least significant bits implies that it is sufficient to 
consider df x where the binary representation of x is 



u n—m—u m—v v terms 


Here Eqs. (22) and (Cl) give 

m 

dfx = 55 q k p m ~ k 55 d J'(*»-*i*-*i 1 -*i) (C3) 

k -0 

with the same limits for the inner sum. These limits can 
be modified by observing that dj( Xn ...x ; ...k !i ...x 1 ) only de¬ 
pends on the number of zeroes and ones in the binary 
representation of x together with the number and ar¬ 
rangement of flips. Suppose that, of the k flips, l occur 
on the rightmost cluster of v bits. Thus k — l occur on 
the second to right cluster of m — v bits. Then the bit 
representation, x n .. . xq .. .xq ... Xi, contains u zeroes 
from the leftmost cluster, 0 from the next cluster to the 
right, [k — l) from the second to right cluster and v — l 
from the rightmost cluster. Thus the total number of 
zeroes in this bit string are u + v + k — 21 and 

dj(x n ...xi k ...xq ...x\) = d u +v+k—2l- (C4) 


Appendix C: Diagonal terms of pf 


Let the binary representation of x be x = x n ... X\. 
Then d{ x = (x|/5f|x) and Eq. (B5) give 


. . Xi l/iprep diag|Xn . . . Xl k . . . Xq . . . Xf) 


(Cl) 


The exact location of the flips within each cluster is irrel¬ 
evant. Thus within the inner sum of Eq. (C3) there are 
(I)Cfc-D wa y g of attaining d u+v+k - 2 i, consisting of (") 
ways of flipping l bits in the rightmost cluster together 
with ( m k Zi) ways of flipping k — l in the second to right 
cluster. Finally not every value of l is possible and these 
are constrained by 0 l ^ v (at most v flips in the right¬ 
most cluster) and 0 ^ k — l^ m — v (at most m — v flips 
in the second to right cluster). Taken together these give 
max (0, v — to + k) ^ l ^ min (k, v). 

Combining these yields Eq. (30). 


Appendix D: Quantum Fisher Information 


We shall show that 
2 


Hx = 


dL ~ 

—2mA 2m-1 df x c 2 


d( x (d t l + m 2 X 2m 2 c 2 (x) 'j 




(Dl) 


and the result will follow from Eq. (31). The expression 
for H x comes from the eigenvalue decomposition of pf x . 
The eigenvalues are of pt x 


p± . df x i A ('j(x) 
and the associated eigenvectors are 

|</>±) := (|z) =F i\N - x))/V2. 


(D2) 


(D3) 


The fact that the eigenvalues are independent of A means 
that the second term on the right hand side of Eq. (8) is 
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zero and the first term then yields 


where 


= 


1 


d{ x + X m Cj^ 

1 


(^df x "b TflX 


+ d,,-\-c 1M f*---*"' 1 '*)) - (D4) 


Straightforward algebra yields the desired result. 

The expression for the quantum Fisher information can 
be rewritten as a sum over u and v. The only relevance 
that each value of x has for the summand is that it de¬ 
termines u, v and j(x) = u + v. Typically many values of 
x will yield the same combination of values of u and v. 
To account for this, note that the sum in Eq. (32) runs 
from x = 00... 0 to x = 01... 1 (in binary notation). 
The fact that the leftmost bit is always zero gives two 
distinct cases. 

In the first, where m < n, this bit structure reveals 
that 1 Sj u ^ n — m and 0 ^ v Sj m. For any choice of u , 
there are f” - ™-/ 1 ) ways in which the zeroes can appear 
in the leftmost n — m bits, given that the leftmost single 
bit must be zero. There are independently (™) ways 
in which the zeroes can appear in the rightmost m bits. 
Thus there are ( n ~l^ -1 ) (™) values of x in the range of the 
sum which result in exactly the same pair of values u and 
v. For each the contribution of the summand is identical. 
This means that the terms within the summand can be 
expressed in terms of u and v (rather than x) and the 
result sum reduces to that of Eq. (33). 

In the second, where m = n, clearly u = 0 and 1 ^ 
v ^ n — 1. A similar counting argument reveals that 
there are ( 1, ~ 1 ) ways in which the zeroes can appear in 
the rightmost n — 1 bits, since the rightmost bit must be 
zero. A similar argument then results in Eq. (34). 


Appendix E: Low polarization correlated state 
protocol quantum Fisher information 


Tn ^max 

E k m-k \ ' 

q p 2^ 

k— 0 1=1 min 

and 





E 


v \ / m — v 


k-l 


k —0 l — l min 

[(2it + 2v + 2k — Al — n) 2 — n\ 


x 


2 "+! 


Both a and /3 are independent of r. Then 

1 


a = 


(E4) 


(E5) 


(E6) 


which can be demonstrated indirectly by noting that, 
when r = 0, each diagonal term of the final density op¬ 
erator equal a. However, when r = 0, po = //2 and thus 
pi = 1/2 n . The channel does not alter this and thus, 
when r = 0, p{ = 1/2”. It follows that a = 1/2”. This 
fact can independently be verified by explicitly evaluat¬ 
ing the three summations in a. 

Note that these imply 



(E7) 


Then Eqs (El), (E3) and (E6), yield that for the terms 
in Eqs. (33) and (34) 


r2 \2m 2 
a i ~ A c u+v 


d{ (df + to 2 A 


2 \ 2m—2 J2 


u u+v 


1 

2^’ 

d{m 2 X~"‘ ~c. 
m 2 A 2r ”- 2 


2 \ 2m—2 „2 

1A+1? 


2 3r 


(2 u + 2v — n) r 2 




d(3 2(2 u + 2v — n) 
d\ 2^+i 


Thus if to < n, Eqs. (33) becomes 


The correlated state protocol quantum Fisher informa¬ 
tion, Eqs. (33) and (34), requires Cj and df (u,v), which 
in turn requires dj. Eqs, (23) and (24) yield, to second 
order in r, 


c, * (2j '- r and 

3 2” 

(El) 

j _ 1 , r2 [( 2 / - n) 2 - n] 

j ~ 2 n 2 n+1 

(E2) 


The latter of these gives 

df ( u , v) = a + f3r 2 (E3) 


m n—m 


iWA)« 2^ E 

0 u— 1 

TO 2 A 2m_ 2 


n — m — 1 
u — 1 


22n+l 


2 3 ” 


(2 u + 2v — n) r 2 


2 m 2 X 2m 2 fn — m — lA 
r 2”- 1 ^ l u — 1 ) 

U= 1 X ' 

m / \ 

x f m J (2 u + 2 v — n). (E8) 

v=o X v / 


The two summations in Eq (E8) give n2” -1 and this 
yields the result of Eq. (39). A similar derivation gives 
the same result when m = n. 
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